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& \ Abstract 



The neutron-star Low-Mass X-ray Binary Aquila X-l was observed seven times in total with the Suzaku 
^vq ■ X-ray observatory from September 28 to October 30 in 2007, in the decaying phase of an outburst. In 

order to constrain the flux-dependent accretion geometry of this source over wider energy bands than 
employed in most of previous works, the present study utilized two out of the seven data sets. The 0.8-31 
keV spectrum on September 28, taken with the XIS and HXD-PIN for an exposure of 13.8 ks, shows an 
■ absorbed 0.8-31 keV flux of 3.6 x 10 -9 erg s _1 cm -2 , together with typical characteristics of the soft state 

of this type of objects. The spectrum was successfully explained by an optically-thick disk emission plus a 
Comptonized blackbody component. Although these results are in general agreement with previous studies, 
the significance of a hard tail recently reported using the same data was inconclusive in our analysis. The 
spectrum acquired on October 9 for an exposure of 19.7 ks was detected over a 0.8-100 keV band with the 
XIS, HXD-PIN, and HXD-GSO, at an absorbed flux of 8.5 x lO" 10 erg s" 1 cm" 2 (in 0.8-100 keV). It shows 
characteristics of the hard state, and was successfully explained by the same two continuum components 
but with rather different parameters including much stronger thermal Comptonization, of which the seed 
photon source was identified with blackbody emission from the neutron-star surface. As a result, the 
accretion flow in the hard state is inferred to take a form of an optically-thick and geometrically-thin disk 
down to a radius of 21 ±4 km from the neutron star, and then turn into an optically-thin nearly-spherical 
, hot flow. 

00 ' Key words: accretion, accretion disks - stars: neutron - X-rays: binaries. 



1. Introduction blackbodies with a range of temperature, and is con- 
sidered to approximate the integrated emission from 

A neutron-star Low-Mass X-ray Binary (NS-LMXB) an optically-thick standard accretion disk (Shakura & 

consists of a weakly magnetized neutron star and a low- Sunyaev 1973). The validity of this model composition 

mass star with the mass typically less than that of the Sun, has been confirmed repeatedly by subsequent studies (?, 

M . NS-LMXBs have been studied from the 1980's and e.g.)]Makishimal989,Asail998,Lin2010,Takahashi2011. 

I were found to dominate the X-ray luminosities of galaxies For comparison, spectra of BHBs in the soft state are 

5^ ■ without strong Active Galactic Nuclei (Makishima et al. even softer (White & Marshall 1984; Makishima et al. 

1989). Like Black Hole Binaries (BHBs), NS-LMXBs are 1986), because black holes lack a solid surface which 

known to exhibit soft and hard spectral states, when their would emit the relatively hard blackbody component, 

mass accretion rate is high and low, respectively (White As represented by figure 1, spectra of NS-LMXBs in 

& Mason 1985; Mitsuda & Tanaka 1986; Mitsuda et al. the hard state are roughly power-law-like with a pho- 

1989). The two types of spectra are visualized in figure ton index of T ~ 2, and often extend up to 100 keV (?, 

1; it compares a soft- and hard-state spectrum of the NS- e.g.)]Tarana2006. Hence in the hard state, broad-band 

LMXB Aquila X-l which we derived by analyzing archival spectroscopy becomes more essential in qualifying major 

RXTE PCA and HEXTE data. emission components. In addition, the luminosity of an 

In the soft state, the spectrum is bright in a soft NS-LMXB in the hard state is not much higher than 

band below 10-15 keV, and falls off steeply towards ~ 10 36 erg s _1 . As a result, the study of NS-LMXBs 

higher energies. It was shown with Tenma that the in the hard state is often limited by the available hard 

soft-state spectra of NS-LMXBs can be expressed by X-ray instrumentation, and has been making a consider- 

the sum of a blackbody component from the NS surface able progress recently thanks to high hard X-ray sensitiv- 

and a disk blackbody component (Mitsuda et al. 1984); ities of BeppoSax, RXTE, INTEGRAL, and Swift. It has 

the latter is defined as a particular superposition of been reported that the spectra can be expressed by the 
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Fig. 1. Background-subtracted vFv Spectra of Aquila X-l 
in the soft state (black) and hard state (red), derived by an- 
alyzing archival RXTE PCA and HEXTE data obtained on 
1998 March 12 and 1999 September 2, respectively. 

sum of an optically thick emission in the soft band and 
a Comptonized emission spanning nearly the whole band 
(Lin et al. 2007; Cocchi et al. 2011; Tarana et al. 2011). 
Also, detailed timing studies are revealing differences be- 
tween NS-LMXBs and BHBs focusing on their hard X-ray 
properties (Titarchuk & Shaposhnikov 2005). However, 
we still need to uniquely associate the spectral emission 
components to actual accretion flows, including the loca- 
tion of the optically-thick emitter, the source supplying 
seed photons of Comptonization, and the geometry of the 
Compton corona. Since Suzaku (Mitsuda et al. 2007) al- 
lows us to acquire the spectra of hard-state LMXBs over a 
very wide (from < 1 keV to > 100 keV) energy band and 
hence to simultaneously quantify the proposed emission 
components, we expect to achieve significant progresses 
on these issues (?, e.g.)]Lin2010. 

In the present research, we focus on the NS-LMXB 
Aquila X-l, and investigate its soft and hard states in 
detail utilizing archived Suzaku data. Aquila X-l is a re- 
current transient with a 2-10 keV luminosity varying by 
a factor of > 200 through its outbursts. It is also an emit- 
ter of Type I X-ray bursts (Koyama et al. 1981), and its 
distance is estimated as 4.4-5.9 kpc by assuming that the 
peak luminosity of its Type I X-ray bursts equals to the 
Eddington luminosity for a 1.4M neutron star (Jonker & 
Nelemans 2004). It was observed by Suzaku seven times 
in 2007, all covering a declining phase of a single outburst. 
Although these data sets were already utilized by Raichur 
et al. 2011, hereafter RMD11, they did not use the HXD- 
GSO data even for the hard state, and hence the available 
information was limited to < 70 keV with relatively large 
errors above ~ 50 keV. We hence try utilizing the GSO 
data for the hard state, taking into account the latest cal- 
ibration results (Yamada et al. 2011). 



2. Observation and Data Processing 

Aquila X-l was observed with Suzaku from September 
28 to October 30 in 2007 for seven times. The log of these 
observations is shown in table 1. Figure 2 is a soft-band 
(1.5-12 keV) light curve obtained with the RXTE ASM. 
Of these 7 observations, RMD11 analyzed the first four. 
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Fig. 2. A soft band (1.5-12 keV) light curve of Aquila X-l 
from 2007 September 1 to 2007 October 31 obtained with the 
RXTE ASM. Arrows indicate the 7 Suzaku observations. 
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Fig. 3. Characterization of the 7 Suzaku observations in 
terms of PIN significance and the signal hardness, the lat- 
ter defined by the HXD-PIN vs XIS count-rate ratio. See 
text (subsection 2.1) for the definition of the two quantities. 
Vertical error bars are statistical only. 



2.1. Characterization of the 7 data sets 

By quickly analyzing the seven Suzaku datasets pro- 
cessed through HEADAS (version 6.9), we obtained the 
source count rates as given in table 1. The processing 
details are described in subsection 2.2. The derived XIS 
count rates for the first four data sets agree with those in 
RMD11, except for a factor 2 difference because they refer 
to a single XIS sensor (though not explicitly stated) while 
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Table 1. Suzaku observations of Aql X-l in 2007. 



No. 


Date* 


ObsID 


ExpT 


Count rate* 








(ks) 


XIS 


HXD-PIN 


1 


9/28 


402053010 


13.8 


275.4±0.4# 


1.12±0.01 


2 


10/03 


402053020 


15.1 


28.12±0.04 


0.74±0.01 


3 


10/09 


402053030 


19.7 


33.18±0.04 


0.89±0.01 


4 


10/15 


402053040 


17.9 


25.85±0.04 


0.73±0.01 


5 


10/19 


402053050 


17.9 


5.67±0.02 


0.13±0.01 


6 


10/24 


402053060 


21.4 


0.221±0.003 


0.041±0.005 


7 


10/30 


402053070 


17.5 


0.172±0.003 


0.056±0.006 



* Date in 2007. 

t Net exposure per XIS sensor. 

* In ct s _1 after subtracting the background. The XIS rate is in 0.8-10 keV, with 
the two FI sensors (XIS0, XIS3) summed. The HXD-PIN rate refers to a 12-50 
keV range, including the CXB contribution of ~ 0.02 ct s _1 . The quoted errors 
are statistical la limits. 

* Corrected for pile-up effects (see text). 



we employ the XIS0 plus XIS3 rates. Our HXD-PIN rate 
is comparable too, but systematically lower by ~ 0.03 ct 
s _1 , particularly in the 2nd-4th observations, than those 
of RMD11: this is presumably because they used a differ- 
ent (non specified) energy band than the 12-50 keV band 
which we used. 

To grasp the source characteristics in the seven obser- 
vations, we calculated, for each of them, spectral hardness 
and HXD-PIN significance. The former is defined as the 
HXD-PIN vs XIS count-rate ratio referring to table 1; the 
latter is given as the background-subtracted 12-50 keV 
PIN signal counts divided by square root of the total PIN 
counts (including the NXB) therein, and provides rough 
measure of the hard X-ray data quality. Figure 3 plots 
the hardness against the PIN significance for the seven 
data sets. The hardness ratio is considerably higher in 
the other observations, of which at least the 2nd, 3rd, and 
the 4th ones were acquired in the hard state according to 
RMD11. Among them, the 3rd observation has the high- 
est PIN significance with hardness comparable to those 
of the 2nd and 4th. Therefore, we have selected for our 
detailed analyses the 1st and 3rd observations, as repre- 
sentative cases of the soft and hard states, respectively. 
The other data sets will be analyzed in a future work, 
including particularly the 6th and 7th ones which are in- 
ferred to exhibit very hard spectra. 

2.2. Data Processing 

2.2.1. XIS data processing 

On September 28 when the source was in the soft state, 
the XIS was operated with 1/4 window and 0.5 s burst op- 
tions. We selected those XIS events of which the GRADE 
was 0, 2, 3, 4, or 6. Time intervals after exits from 
the South Atlantic Anomaly were set to be longer than 
436 s. Since the XIS count rate was so high as to need 
pile- up corrections, we discarded the image center (< V) 
when accumulating on-source events as described above. 
Background events were obtained from a larger annulus, 
which spans a radius range of 4' to 6'. A pile- up corrected 
and background-subtracted XIS FI (XIS0 + XIS3) count 



rate was 137.7 ct s . The count rate in table 1 was ob- 
tained by restoring the counts to be observed within the 
excluded central region (< V) with the help of the point- 
spread function (Serlemitsos et al. 2007); this process is 
expected to be accurate to ~ 15% (Yamada et al. 2011). 
Due to the burst option, the net exposure of XIS0 and 
XIS3 was 3.5 ks each, which is a quarter of 13.8 ks shown 
in table 1. 

In the other observations including the October 9 one, 
when the object was presumably in the hard state, the XIS 
was operated with the 1/4 window option, but without 
incorporating the burst option since the expected count 
rate was not very high. We accumulated the on-source 
events within a circle of 3' radius, because the pile-up 
correction was not required. The background events were 
taken from an annulus with the inner and outer radii of 
3' and 5', respectively. 

Light curves of the processed XIS data of the two obser- 
vations (the 1st and 3rd), employed for our detailed study, 
are shown in figure 4, panel (a) and (c). In either observa- 
tion, no burst-like events are seen and the time variation 
is less than ~ 10%. Therefore, for our spectral analysis, 
we utilized time- averaged data, accumulated over a 0.8-10 
keV range. To avoid the instrumental silicon K-edge and 
gold M-edge, where calibration uncertainties are large, we 
excluded two energy ranges, 1.7-1.9 keV and 2.2-2.4 keV, 
respectively. 

2.2.2. HXD data processing 

In all observations, the HXD data were screened by con- 
ditions that the time intervals after exits from the South 
Atlantic Anomaly should be longer than 500 s, the eleva- 
tion of the target above the Earth limb should be > 5°, 
and the geomagnetic cutoff rigidity should be higher than 
8 GV. The Non X-ray Background (NXB) spectra were 
created from a fake-event file provided by the HXD team, 
and then subtracted from the on-source spectra. In any 
observation, the Cosmic X-ray Background (CXB) was 
not subtracted from the data at this stage but included 
as a fixed model in the following spectral fits (section 3). 
The HXD-PIN count rates in table 1 were calculated be- 
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fore removing the CXB contribution. 

Analyzing the soft-state (September 28) data, we de- 
tected positive source signals (after subtracting the NXB 
and CXB) in the energy band of 12-50 keV, but we quote 
the source detection conservatively in 12-31 keV, where 
the signal exceeded 5.6% of the NXB which is a typi- 
cal 2.5a error for an exposure of ~ 14 ks (?, calculated 
from)]Fukazawa2009. Further examination on this point 
is given in subsection 3.1.4. The HXD-GSO data were 
consistent, within statistical errors, with the background 
over the whole energy range (50-600 keV). In the October 
9 observation, which we have selected as the representa- 
tive of the hard state, the HXD-PIN and HXD-GSO sig- 
nals were successfully detected over 12-50 keV and 50-100 
keV, respectively. The 50-100 keV HXD-GSO count rate 
was 0.19 ±0.03 ct s _1 . This is significant, even consid- 
ering systematic NXB error which is ~ 0.6% (la) for the 
HXD-GSO case, or ~ 0.045 ct s _1 for an exposure of ~ 10 
ks (Fukazawa et al. 2009). HXD-PIN light curves of the 
selected two observations are shown in figure 4 panel (b) 
and (d). 

3. Spectral analysis 

3.1. The soft state 

In performing quantitative model fitting to the 
September 28 spectra using XSPEC (version 12.6.0), we 
compensated for the reduced data accumulation region 
(discarding < 1') by properly calculating an XIS arf that 
excludes the central circle within V . To fit the HXD-PIN 
spectrum, the CXB was included in the model as a fixed 
component, expressed as a function of energy E as 

CXB(^) = 9.41xlO-(^)"exp(- iU | v ) (1) 

where the unit is photons cm _2 s _1 keV _1 FOV _1 . 
Although the CXB surface brightness in the HXD- 
PIN range has an uncertainty of about ±15% (?, 
e.g.)]Turler2010, this is by an order of magnitude smaller 
than typical statistical errors in our HXD-PIN spectra, 
and is hence negligible. The cross- normalization difference 
between the XIS and the HXD can be considered by mul- 
tiplying the model for HXD-PIN with a constant factor of 
1.158, which was established accurately through calibra- 
tions using the Crab Nebula data obtained on many differ- 
ent occasions (Ishida et al. 2006; Ishida et al. 2007; Maeda 
et al. 2008). Its uncertainty, ±0.013, is smaller than the 
statistical errors associated with the present HXD-PIN 
data, but later we allow it to vary. 
3.1 .1. Fit with diskbb + blackbody 

Figure 5 shows the XIS (FI) and HXD-PIN spectra ob- 
tained on September 28, presented after subtracting the 
NXB, but without removing the instrumental responses. 
They indeed bear typical features of the soft state. We 
jointly fitted them with a canonical model (Mitsuda et al. 
1984; Mitsuda et al. 1989; Makishima et al. 1989), which 
consists of a soft disk blackbody component and a hard 
blackbody component. Hereafter, the CXB model as de- 
scribed above is included as a fixed component for the 



HXD-PIN spectra. The spectra were approximately re- 
produced as shown in figure 5(a) and (b), in terms of 
the blackbody temperature of Tbb ~ 1.9 keV, its radius of 
i?bb ~ 2 km, the inner disk radius of R[ n ~ 11 km, and the 
temperature therein T- in ~ 0.8 keV. However, the fit was 
not statistically acceptable with xi — 3.24(146). Since xi 
is larger than 2.0, errors cannot be assigned to the pa- 
rameters, and the residuals are presented in the form of 
data vs model ratio rather than in units of %. This model 
failure results from noticeable negative residuals around 
7.0-9.0 keV and positive ones above 20 keV. 
3.1.2. Fit with an optically -thick component plus its 
Comptonization 

The positive fit residuals seen in figure 5(a) above ~ 15 
keV are suggestive of inverse- Compton effects, working 
possibly on blackbody as reported in previous studies 
(Mitsuda et al. 1989; Gierlihski & Done 2002; Lin et al. 
2007; Cocchi et al. 2011; Tarana et al. 2011; Takahashi et 
al. 2011). We hence tried, as a first-cut attempt, mod- 
els involving a single optically-thick emission, assuming 
that its fraction is Comptonized while the rest is directly 
visible. As the thermal Comptonization model, we em- 
ployed nthcomp (Zdziarski et al. 1996; Zycki et al. 1999) 
here, because it can select between a blackbody and a 
disk blackbody as the seed photon source, and can han- 
dle a relatively wide range of optical depths. Thus, two 
models were constructed, diskBB+nthcomp(diskBB) and 
bbody+nthcomp(bbody) , where the seed photon source to 
nthcomp is indicated in parentheses. The photon temper- 
ature was left free but constrained to be the same between 
the direct emission and the Comptonized one, while their 
normalizations were both left free, as well as the optical 
depth of nthcomp. As a result, the former model gave 
T in = 0.59±g;g| keV with the direct fraction of 0.6 ±0.1 
(60% being direct), while the latter model T bb = 0.381^2 
keV with the direct fraction of 0.3 ±0.1. However, as pre- 
sented in figure 5(b) and figure 5(c), neither fits were suc- 
cessful yet since they gave xt{v) = 1.56(146) and xt( u ) = 
1.60(146), respectively. These fits were not improved at 
all even when the cross-normalization between the XIS 
and HXD-PIN was allowed to vary; xt(y) = 1.57(145) and 
X 2 » = 1.61(145). 

The above fit failures are due, at least partially, to 
the noticeable negative residuals seen in figure 5 (pan- 
els b and c) at energy ranges of 7.0-7.5 keV, 8.5-9.0 
keV, and possibly 1.5-1.6 keV. To account for the two 
higher-energy features, we hence added a neutral Fe K- 
edge and an ionized Fe K-edge, with their energies both 
left free. To eliminate that around 1.5-1.6 keV, which 
could arise from calibration errors of the Al K-edge in 
the X-Ray Telescope (Serlemitsos et al. 2007) of Suzaku, 
another edge absorption was added with its energy fixed 
to 1.56 keV. Then, the diskBB+nthcomp(diskBB) and 
bbody+nthcomp(bbody) fits were significantly improved 
to xliy) = 1-34(141) and 1.27(141), respectively. Fit resid- 
uals of the former case are presented in figure 5(d). Taking 
the latter case for example, the two Fe-K edge energies 
were obtained at 7.2 ±0.1 and 8.8 ±0.2 keV, correspond- 
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Bin time: 120.0 s 




15 20 25 

Time (ks) 

Start Time 14371 15:41:47:184 Stop Time 14372 1:21:47:184 



20 
Time (ks) 

Start Time 14382 5:29:09:184 Stop Time 14382 16:31:09:184 



Fig. 4. Light curves of Aquila X-l with a 120 s binning obtained on 2007 September 28, in 0.8-10 keV with the XIS (XIS0 plus 
XIS3: panel a) and 12-50 keV with HXD-PIN (panel b). Panels (c) and (d) are the same as (a) and (b), respectively, but obtained 
on 2007 October 9. The HXD-PIN light curves include the CXB, but not the NXB. The rate in panel (a) is ~ 50% of that in table 
1, due to the exclusion of the events within V '. 



ing to those of neutral Fe and He-like or somewhat lower 
ionization species, with the associated optical depths of 
0.10 ±0.04 and 0.14 ±0.06, respectively. Thus, the edges 
were confirmed to be significant and their parameters are 
generally reasonable, although neither model has achieved 
adequate fit goodness yet. 

3.1.3. Fit with a disk blackbody and a Comptonized 
blackbody 

Since the models with single seed photon population 
(panels b, c, and d of figure 5) were found to be unfavor- 
able, we proceed to double seed photon population mod- 
els, in which the direct emission and the Comptonized 
emission have different optically-thick radiation sources. 
As the most natural extension from the classical model 
utilized in subsection 3.1.1 and presented in figure 5(a), 
a model of the form diskBB+nthcomp(bbody) was con- 
structed, in which the two optically-thick components 
were allowed to take different temperatures, T{ n and Tbb- 
This describes a configuration where the disk emission is 
directly visible while the entire blackbody emission from 
the NS surface is Comptonized. When the three edges are 
not considered, the fit gave xt{v) = 1.44(145). Although 
this is not yet acceptable, it is significantly better (by 
Ax 2 = -18.4 to -24.1 for Av = -1) than those obtained 
in the previous subsection under the single-source model- 
ings. 

When the three edges are incorporated, the fit has be- 
come fully acceptable to xt{v) = 1.15(140), with a null- 
hypothesis probability of 0.1. The obtained fit is shown 



in figure 5 (panels e and f), and the model parame- 
ters are summarized in table 2; the results virtually re- 
mained the same even when leaving free the XIS vs. HXD 
cross-normalization. The nthcomp model specifies the 
Comptonization via two parameters, namely, the coronal 
electron temperature T e , and the spectral slope T c . In ta- 
ble 2, we have converted this T c (obtained as 3^) to the 
optical depth of r ~ 6, using relations as 



(2) 




(Lightman & Zdziarski 1987). For reference, the Compton 
^/-parameter, y= (4/cT e /ra e c 2 ) x max(r,r 2 ), becomes 0.8 ± 
0.1. 

In table 2 and hereafter, the innermost disk radius R[ n 
has been derived from the diskBB model, assuming that 
its radius parameter is in fact equals to £ft 2 i?i n Vcosz. 
Here, £ = 0.412 is a correction factor for the inner bound- 
ary condition (Kubota et al. 1998; Makishima et al. 2000), 
n = 1.7 is the standard color hardening factor (Shimura 
& Takahara 1995), and i is the disk inclination which we 
assumed to be 45° for simplicity. Thus, the estimated 
emission-region radius can be reasonably identified with a 
fraction of the NS surface. 

Just for reference, we examined another possible model 
configuration, namely bbody + nthcomp (diskBB) . Like 
in the case of table 2, all model parameters (again ex- 
cept edgel energy at 1.56 keV) were left free. It yielded 
xl(v) = 1.53(145) and 1.30(140), without and with the 
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Fig. 5. (a) Simultaneous fitting of the XIS and HXD-PIN spectra in the soft state and its data-to-model ratio, using a diskBB (red) 
and a bbody (blue) model, (b) Residuals of the fit using a single diskBB plus its Comptonization. (c) The same as (b) but the 
diskBB source is replaced with a blackbody. (d) The residuals incorporating three absorption edges to (b). (e) A fit with diskBB 
plus nthcomp (bbody), and its residuals, (f) The vFv form corresponding to panel (e). The error bars represent statistical la level. 



edges, respectively. The fit was not improved [xli^) = 
1.29(139)] even when the cross- normalization was left free. 
We no longer consider this alternative modeling, since 
these fits are significantly worse than those obtained with 
diskBB+nthcomp (bbody) , and are not acceptable. 
3.1.4. Significance of the hard tail 

Analyzing the soft-state Suzaku data on September 28, 
RMD11 reported the detection of a hard-tail component 
with a very flat slope of T ~ 0, which appeared in the 
30-70 keV range of the HXD-PIN spectrum. To ex- 
amine its reality, we extended the upper energy bound 
of the HXD-PIN data from the previously employed 31 
keV to 70 keV, by adding five data points. As shown 
in figure 6 (a), these additional data points indeed ex- 
ceed the best-fit diskBB+nthcomp (bbody) model (with 
the three edges), even when the model is refitted. Then, 
by adding a PL with T = (fixed), the excess was ex- 
plained away as seen in figure 6 (b), and the fit good- 
ness changed slightly from xt(i/) = 1.15(145) to 1.12(144) 
with Ax 2 = -5.9 for Ai/ = -1. Within 90% errors, the 
obtained PL normalization, (9 ± 6) x 10 -6 , agrees with 
(2.1 ± 1.5) x 10" 5 reported by RMD11. According to an 
F-test, the probability for this fit improvement to arise 
by chance is about 2%, even though the statistical devi- 



ation of the additional five data points from the best-fit 
diskBB+nthcomp(bbody) model yields a relatively low re- 
duced chi-square of xt — l-25(^ = 5). When allowed to 
vary, T was not well constrained, and floated over a rather 
large range, e.g., 0-3.0. In any case, the other fit param- 
eters were virtually unaffected, and remained consistent 
with table 2. 

The above estimates considered only statistical errors. 
However, we must also consider systematic uncertainties 
in the HXD-PIN data points. The NXB modeling in the 
40-70 keV range of HXD-PIN is known to have la sys- 
tematic errors of 2.8% and 1.8%, for an exposure of 10 
ks and 20 ks, respectively (Fukazawa et al. 2009). Since 
the present data have an exposure of 14 ks, we may sim- 
ply estimate the la NXB error as 2.2%, as the geometric 
mean of 2.8% and 1.8%. Similarly, we assign a la system- 
atic NXB error of 1.2% below 40 keV. Figure 6(c) shows 
the same spectrum as figure 6(b), but in the vFv form, 
in comparison with these la NXB uncertainties. When 
this systematic error is included in the fit goodness evalu- 
ation, the F-test chance probability for the additional PL 
component increased to ~ 7%. Therefore, the hard tail is 
significant at the 90% confidence level, but not at 95%. 




Fig. 6. The same soft-state spectra as figure 5, but further including five additional HXD-PIN data points in the 31-70 keV range, 
(a) The diskBB+nthcomp(bbody) fit, which is nearly identical to that in figure 5 (e), but the parameters were re-optimized, (b) The 
same as panel (a), but adding a PL (green dotted line) with a fixed slope of F = 0.0. (c) The vFv form of panel (b), where systematic 
lcr uncertainties of the NXB (see text) are superposed in purple. Again, the error bars of the data points represent statistical la 
range. 



3.2. The hard state 

The October 9 spectra, detected over the 0.8-100 keV 
band with the XIS, HXD-PIN, and HXD-GSO (subsec- 
tion 2.2.2), were analyzed in the same way. As presented 
in figure 3 and reported by RMD11, the spectra on this 
occasion show typical characteristics of the hard state; 
we therefore consider that the source was actually in the 
hard state on this occasion. Since the interstellar absorp- 
tion, corresponding to wabs in table 2, was obtained as 
TVh = 0.36 x 10 22 cm -2 in the soft state analyses where 
the soft band flux is higher, we fixed iVn to this value 
in the hard state analyses. Neither of the three edges 
were incorporated, because the soft band flux was rela- 
tively low in this observation. The model normalization 
was constrained to be the same between HXD-PIN and 
HXD-GSO (Kokubun et al. 2007). 
3.2.1. Fit with a single Comptonized component 

Since hard-state spectra are generally considered to be 
more dominated by Comptonization than the soft-state 
ones analyzed in subsection 3.1, we begin with fitting 



the October 9 spectra with a single Comptonized emis- 
sion model, represented by compPS (Poutanen & Svensson 
1996). The use of this model, instead of nthcomp used 
in subsection 3.1, is because it can more properly han- 
dle relativistic effects (e.g. using the Klein-Nishina for- 
mula instead of the Thomson cross section), and hence is 
more suited to the hard-state spectrum which extends to 
~ 100 keV. Actually, this model was successfully employed 
by Takahashi et al. 2008 and Makishima et al. 2008, in 
the study of GRO 1655—40 and Cygnus X-l, respectively, 
both observed with Suzaku in the hard state. Hereafter, 
a relatively narrow iron line is added to the fitting model. 
As to seed photons of the inverse Compton emission, we 
assumed either a blackbody with free Tbb or a diskBB 
with free T[ n . As shown in figure 7(a), the spectra were 
roughly reproduced by a strong Comptonization operat- 
ing on either bbody or diskBB (though only the former 
is shown in the figure), with Tbb ^ 0.5 or T in ~ 0.6 keV, 
respectively. However both of them were unacceptable, 
mainly due to positive residuals below ~4 keV. 
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Table 2. Parameters of the diskBB plus nthcomp fit for the 

soft state.* 



Component 


Paramater 


Value 


wabs 


N u (lO^cm- 2 ) 


0.36 ±0.01 


edgel 


&edge (KeV J 

optical depth 


i.oo (nxeaj 
0.03 ± 0.01 


edge2 


E edge (keV) 


7.1 ±0.1 




optical depth 


0.14 ±0.04 


edge3 


£7 i (keV) 


8.7 ± 0.2 




optical depth 


0.13 ± 0.06 


diskBB 


Ti„ (keV) 
Rin (km)t 


n 7 o+0.U2 
u - ' °-0.03 

15 ±2 


nthcomp 


T bb (keV) 
T e (keV) 
optical depth 
i?bb (km) 


1 4+^ 

o+28 

6±4 

3.ot° ; 8 7 


Fit goodness 


x» 


1.15 (140) 



* Errors represent 90% confidence limits, 
t After applying the corrections described in the text, and 
assuming a distance of 5.2 kpc. 

When leaving Nn free, the compPS(bbody) fit was 
much improved from y? v {y) = 8.64(186) to 1.85(185), but 
the result was not yet acceptable because of the struc- 
tures below rsj 4 keV. By this additional degree of free- 
dom, the compPS (diskBB) fit improved from xl(v) = 
2.52(186) to 1.16(185). However, the derived value of 
Ah = 0.277to;? 5 x 10 22 cm -2 is not consistent with that 
of the soft state (table 2), and the obtained value of 
Ri n = 4 ± 1 km is too small. Leaving free the XIS vs HXD 
cross normalization did not improve either fit. Therefore, 
the spectrum cannot be reproduced by a single thermal 
Comptonization model, regardless of the choice of the seed 
photon source. More specifically, the data are suggestive 
of the presence of an additional component at the softest 
spectral end, and a reflection hump at ~ 30 keV. 
3.2.2. Fit with an optically-thick component plus its 
Comptonization 

To account for the soft excess left by the single 
compPS fits in figure 7(a), we modified the model, 
like in the soft-state analysis, and considered a case 
where a fraction of the seed photon source (either 
bbody or diskBB) is directly visible. Hence the model 
is bbody±compPS(bbody) or diskBB±compPS(diskBB). 
Again, the direct and Comptonized photons were assumed 
to have the same seed-photon temperature. However, in 
both cases, the fraction of directly seen component be- 
came so small (~0.04 and ~0.001 for diskBB and bbody, 
respectively), that the results were nearly the same as in 
the previous subsection. This is mainly because the soft 
excess, appearing in < 1 keV, is rather softer than the 
Compton seed photon source, of which the photon tem- 
perature is specified as ~ 0.5 keV by the ~ 1 to ~ 4 keV 
part of the spectrum. In short, the spectra cannot be 
reproduced by a single photon source, even allowing its 
fraction to be directly visible. 



3.2.3. Fit with an optically -thick component plus a sep- 
arate Comptonized emission 

After previous works (Christian & Swank 1997; Church 
& Balucihska-Church 2001; Gierlihski & Done 2002), as 
well as the results of subsection 3.2.2, we proceed to a 
fit in which the directly-visible soft thermal emission is 
different from the seed source for compPS. 

As one of such model forms, we tried 
diskBB±compPS(bbody) model, which has essentially the 
same form as our final model for the soft state (except 
Fe-K edges and lines). As shown in figure 7(b), the 
soft-band data below ~ 2 keV have been well reproduced 
(hence with little merit of leaving Ah free), but the fit 
was not yet acceptable with xt{v) = 1.43(184), due to a 
step-like feature in 7-10 keV and a hump at ~ 30 keV. 
These features were not reduced significantly by changing 
the XIS vs HXD cross normalization. Since they are 
suggestive of the presence of a reflection component like 
in the soft-state spectra, we incorporated a reflection 
embedded in compPS. Then, as shown in figure 7(c) and 
table 3, the fit has become acceptable even when Ah is 
fixed. The obtained parameters, including the reflection 
solid angle, are reasonable for the modeling. 

As an alternative attempt, we fitted the spectra with 
the other model form, namely, bbody±compPS(diskBB). 
Like the previous model, this fit was unacceptable with- 
out considering reflection, xt{v) = 2.05(184), even when 
leaving TVh and the cross normalization free. When re- 
flection is incorporated, the fit has become acceptable 
(with Ah fixed), as shown in figure 7 (d) and table 3. 
The final fit goodness is somewhat worse than that of 
the diskBB±compPS (bbody) modeling with the fixed TVh- 
Further comparison between the successful two models is 
presented in section 4.3. 

4. Discussion 

4-1. The two spectral states 

We analyzed two out of the 7 Suzaku data sets of 
Aquila X-l, acquired in an outburst over 2007 September 
to October. On September 28, when the source was in 
the typical soft state as expected from figure 3, the sig- 
nals were successfully detected over 0.8-31 keV with the 
XIS and HXD-PIN. On October 9, in the typical hard 
state (figure 3), the spectrum became much harder, and 
the signal detection was achieved over a broader range of 
0.8-100 keV, with the XIS, HXD-PIN, and HXD-GSO. 
The soft-state luminosity in 2-20 keV and that of the 
hard state in 2-100 keV were 1.5 x 10 3T erg s _1 and 
2.5 x 10 36 erg s _1 , respectively, assuming a distance of 
5.2 kpc and isotropic emission. These values correspond 
to ~8% and ~1% of the Eddington luminosity, respec- 
tively, for a neutron-star mass of 1.4 M . Therefore, the 
present results are consistent with the understanding that 
NS-LMXBs, in a broadly similar way to BHBs, make tran- 
sitions between the soft and hard states at a threshold 
luminosity of several percent of the Eddington luminosity 
(Bloser et al. 2000; Maccarone 2003). Below, we discuss 
the two states separately, referring to the results obtained 
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Fig. 7. Simultaneous fitting of the XIS, HXD-PIN, and HXD-GSO spectra in the hard state, (a) A fit with compPS (bbody) and its 
data-to-model ratio, (b) A fit with diskBB+compPS (bbody) and its ratio, (c) A fit with the same model as (b), but incorporating 
reflection. The middle panel shows residuals, and the bottom one vFv form of the fit. (d) A fit obtained by exchanging the seed 
photon and the optically-thick emission in (c), namely, with the bbody + compPS (diskBB) model. 
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Table 3. Fit parameters of an optically-thick component plus a Comptonized emission for the hard state. 



Component diskBB±compPS(bbody) bbody±compPS(diskBB) 



Nh (10 22 cm- 2 ) 




0.36 (fixed) 


4^+0.05 
U ' 4O -0.07 


0.36 (fixed) 


0.50 ±0.06 


Opt. thick 


model 


diskBB 


diskBB 


BB 


BB 




T bb /T in (keV) 
Rbb/Rin (km) 


0.278 ±0.02 
21±4 


0.23™ 

4-1+I8 


0.16 ±0.01 
41 ±8 


m ^ ^ ,4+0 009 

0.144™? 
100 ±30 


compPS 


seed 


BB 


BB 


diskBB 


diskBB 




T in /T bb (keV) 
Rin/Rbb (km) 
T e (keV) 
r 

refl* (fi/27r) 


0.51 ±0.02 
10 ±2 

35±S 
> 2.5+ 

o.6t° ; 2 


0-48±g;g| 
12 ±3 
37±t 

9 r+0.3 
Z '°-0.2 

7 +0 - 2 


~ — . ^ -lc\ nn^ 

u. ' z ±y_ 004 
4.2 ±0.6 

52.418;5 

-I 70+O.OI 
0.60 ±0.05 


0.68 ±0.04 
4.3 ±0.8 

5515 

i-e±8:S 
o.6t° ;? 


gaussian 


E (keV) 
a (keV) 
EW* (eV) 


6.4 (fixed) 
0.1 (fixed) 

oo+13 
zz -16 


6.4 (fixed) 
0.1 (fixed) 

21 +14 
Zi -16 


6.4 (fixed) 

0.1 (fixed) 
90+15 

zo -16 


6.4 (fixed) 
0.1 (fixed) 
22 ±17 


fit goodness 


xiiy) 


1.02(183) 


1.00(182) 


1.10(183) 


1.04(182) 



* Reflection solid angle. 

t The upper limit is larger than 3, which is the model limit in compPS. 
$ Equivalent width. 



with Nn = 0.36 x 10 22 cm 2 and the cross normalization 
of 1.158 both fixed. 

4.2. Geometry in the Soft state 

Using the parameters derived from the successful fit 
with diskBB±nthcomp(bbody) shown in figure 5(e), we 
obtained a blackbody radius of i? bb = 3.0^0.8 ^ m anc ^ an 
inner disk radius of R m = 15 ±2 km (table 2). According to 
these results, the optically-thick (and geometrically-thin) 
disk is inferred to continue down to close vicinity of the 
neutron-star surface, and the matter then accretes onto 
the surface mainly along an equatorial zone to emit the 
blackbody radiation. The blackbody emission from the 
surface is moderately Comptonized, with a ^-parameter 
of ~ 0.8 (subsection 3.2.3), by a sort of atmosphere or 
corona, where the accretion flow would be rather radia- 
tively ineffective and the internal energy would become 
higher therein. Like in BHBs, the coronal electron tem- 
perature, T e , is considered to be determined via a bal- 
ance between heating by ions and Compton cooling by 
the bbody photons. 

To re-examine the very flat hard-tail component 
claimed by RMD11 using the same data set, we fur- 
ther included five data points in the 31-70 keV range 
(subsection 3.1.4). They indeed exhibit excess above the 
diskbb±nthcomp (bbody) fit (figure 6), and this effect can 
be accounted for by a hard PL (with r ~ fixed) with nor- 
malization consistent with what was reported by RMD11. 
However, the PL component has a rather marginal signif- 
icance (confidence level of 93%), when properly consider- 
ing systematic errors associated with the HXD-PIN NXB 
subtraction. In addition, the suggested phenomenon, with 
r ~ 0, is rather unusual. For these reasons, we refrain from 
further discussion on this issue. 

According to the virial theorem, matter in a standard 
accretion disk must radiate half its gravitational energy 
which is released by the time it reaches R[ n . Assuming 



that the corona is localized in a region around the neu- 
tron star, the remaining half energy, plus the energy 
released from R[ n to the neutron-star radius i?NS> will 
power the blackbody from the neutron-star surface and 
the corona. Therefore, when spherically integrated, the 
disk luminosity is expected to be comparable to that of the 
Comptonized blackbody, assuming R{ n ~ Rns- However, 
the observed fluxes of diskBB and nthcomp (bbody), de- 
noted Fdisk and F bb respectively, will depend on our 
viewing angle in different ways. As noticed above, the 
blackbody-radiating area would have a shape of a thin and 
short cylinder along the equatorial region of the neutron 
star, and so would be the corona. Then we can estimate 
the inclination angle 0\ of the disk as 



(Mitsuda et al. 1984). The observed ratio of F bb (in 0.1- 
100 keV) to F disk (in 0.1-50 keV) is 1 : 2.3, which con- 
strains the inclination angle as ~ 34°. This value is not 
so different from the value of 45°, which we tentatively 
assumed in the data analysis. 

The successful fit incorporated two absorption edges 
that are considered of celestial origin. Their energies 
turned out to be 7.lto'i keV and 8.5 ± 0.4 keV, corre- 
sponding to neutral Fe-K edge and that of highly ionized 
Fe-K, respectively. However, the data lack associated Fe- 
K fluorescence lines (with an upper-limit equivalent width 
of 10 eV assuming them to be narrow), like in a previous 
study (D'Ai et al. 2006); the absence of lines might be 
attributed to the geometry of the absorber. 

4.3. Geometry in the Hard state 

The hard-state spectra have been explained success- 
fully by either a diskBB± compPS (bbody) model or a 
bbody± compPS (diskBB) model. Based on the fit good- 
ness alone, we cannot tell which of the two model forms is 
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more appropriate (subsection 3.2.3). Let us then examine 
the two solutions for their physical appropriateness. In 
the diskBB+compPS(bbody) modeling, the derived radii 
are = 10 ± 2 km and R[ n = 21 ± 4 km. The former is 
nearly the same as that of a neutron star, and the latter 
is larger than that; this is physically reasonable. In the 
bbody+compPS(diskBB) modeling, in contrast, the radii 
are = 41 ± 8 km and R[ n = 4.2 ± 0.6 km, which are 
not physical since these would mean that the inner edge 
of the disk is smaller than the neutron star radius. Thus, 
the diskBB+compPS(bbody) model, which has the same 
form as our final soft-state model, is clearly more favored. 

The same October 9 data, together with the other two 
hard-state data sets (2 and 4 in figure 3), were already 
analyzed by RMDll. They used several models which are 
qualitatively similar to ours, each consisting of a softer 
optically-thick emission and a harder Comptonized com- 
ponent. Specifically, their Mia model has a form of 
diskBB+nthcomp(diskBB), and their Mlb model is de- 
scribed as diskBB+nthcomp (bbody). Although their Mlb 
model is essentially the same as our final solution, RMDll 
concluded Mia to be the best model. Furthermore, pre- 
sumably due to differences in the employed Compton 
codes and in the Fe-K line modeling, our results and theirs 
show some discrepancies as listed below. 

1. Our successful final fit with diskBB+compPS (bbody) 

to the hard-state data has given a reasonable in- 
ner disk radius as R[ n = 21 ±4 km (table 3), while 
RMDll found, with their favourite Mia model and 
assuming the same distance, too small a value of 
R in = 6.7 ±0.2 km. 

2. The value of T- in ~ 0.3 keV we measured in the hard 
state is lower than that (~ 0.7 keV) in the soft state, 
while RMDll reported (again with the Mia model) 
opposite behavior, which would be inconsistent with 
the change in the mass accretion rate. 

3. The values of TVh which we obtained are consistent 
with being the same between the two states, while 
this does not apply to the Mia results of RMDll. 

4. Even if we compare our results with the Mlb fit by 
RMDll (which has the same model form as ours but 
is not favored by RMDll), they obtained too large 
a value of i?bb = 18 d= 5 km, while our estimate of 
i?bb = 10 =b 2 km (table 3) is fully consistent with 
the neutron-star radius. Furthermore, we found in 
both states T- in < Tbt>, while this is not necessarily 
the case in the Mlb results of RMDll as they admit. 
(No information as to these points is available with 
their Mia modeling.) 

Considering these, as well as the appropriateness of 
compPS to very hard spectra, we consider that our hard- 
state results provide a better description of the October 9 
data. 

Although our final spectral model for the hard state 
has the same composition as that in the soft state (ex- 
cept the edges and line), the model parameters are con- 
siderably different, as we see below. The standard disk 
is truncated at a larger radius (~20 km) than in the 



soft state. Inside that radius, the flow becomes optically 
thin and geometrically thick. This optically-thin disk it- 
self would be identified with the Compton corona, with 
a much larger ^-parameter of ~ 2 than in the soft state 
(y ~ 0.8). Unlike that of the soft state, the accretion flow 
around the neutron star is considered to be rather spher- 
ical, because nearly the whole neutron-star surface is in- 
ferred to be emitting the blackbody seed photons, which 
are then inversely- Comptonized by the corona. Although 
it is still possible that the standard disk is partly covered 
with the corona and part of its emission is Comptonized, 
the Comptonization of the blackbody should be rather 
dominant because the reflection solid angle is not so large 
(ft/27r ~ 0.6). 

The observed flux ratio of compPS (bbody) to diskBB 
is about 6:1. Thus, the relative dominance of the two 
components has reversed as compared to the soft state. 
This is explained, at least qualitatively, by the following 
two effects. One is the increase in i? in , which would reduce 
^disk (as oc i?^ 1 if the total luminosity is conserved) and 
enhance F^. The other is the geometrical change in the 
corona and the blackbody, from the short cylinder in the 
soft state to the nearly spherical configuration in the hard 
state. 

The fit required a narrow line at 6.4 keV, with an equiv- 
alent width of ~ 22 eV. This is interpreted as (nearly) 
neutral Fe-K line, produced when the cool accretion disk 
is irradiated by the very hard continuum. 

In the picture as derived so far, the motion of the ac- 
creting matter is expected to change from Keplerian to 
nearly free-fall like, at a radius of R[ n ~ R c , where R c is 
a typical radius of the corona. Therefore, the infall speed 
of the matter in the corona may be written as 



v(r) = g 



2GM N 



(0 <<?<!), 



(4) 



where r is the distance from the neutron star, G is the 
gravitational constant, M NS is the neutron star mass, and 
g is a numerical factor of order unity. Meanwhile, the 
observed total luminosity is represented as 



-^obs — / 



GM NS M 

Rns 



(0</<l), 



(5) 



where M is the mass accretion rate, and the numerical 
factor / becomes less than 1 when there are outflows or 
when a fraction of the released energy is stored as internal 
energies of the neutron star. This M is also related with 
the radius and electron density n e (r) as 



M = 4:7rr 2 v(r)m p n e (r), 



(6) 



where m p is the proton mass. Using the three equations, 
and performing integration from r = i? NS (=10 km) to r = 
i? c , the optical depth of the corona is represented as 



= cr T h / 



n e (r)dr 



:0.04x <^ 1 



V 10km / 



1 

Tg' 



(7) 
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where a Th is the Thomson cross section. So, if we assume 
Rin ^ Rc ^ 20 km and the product takes a value of 
~ 0.01, we can explain the observed value of r. 

5. Summary 

Utilizing the archival Suzaku data of Aquila X-l in an 
outburst, we successfully obtained high-quality spectra on 
two occasions, and performed quantitative model fitting 
to both spectra. The obtained results are summarized as 
follows. 

1. The 0.8-31 keV spectra obtained on 2007 September 
28 exhibited typical characteristics of the soft state, 
with the 2-20 keV luminosity being ~ 8% of the 
Eddington limit. In contrast, the 0.8-100 keV spec- 
tra on 2007 October 9, much harder as usually seen 
in the hard state, had a 2-100 keV luminosity which 
is ~ 1% of the Eddington limit. 

2. The results of the soft state analyses are generally 
consistent with those from previous studies, includ- 
ing RMD11. There are two emission components; a 
disk blackbody emission from a standard disk, and 
a weakly Comptonized blackbody emission arising 
from the neutron-star surface. 

3. The flat hard-tail component, claimed by RMD11 
based on the same Suzaku observation on September 
28, becomes rather inconclusive when considering 
the systematic background errors. 

4. The hard-state spectra can be expressed by the sum 
of a disk blackbody and a Comptonized blackbody. 
The Compton seed photons are provided by a black- 
body emission from the neutron star surface. 

5. In the hard state, a standard accretion disk is in- 
ferred to be truncated at ~ 20 km, at which point 
the accretion flow turns into an optically-thin nearly 
spherical flow. The blackbody emission from the 
neutron-star surface is strongly Comptonized with 
a ^-parameter of ~ 2. For these two reasons, the 
overall spectrum is much harder than that in the 
soft state, even though the basic emission compo- 
nents are the same. 
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